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Abstract 

We introduce a new method from number fields and codes to construct dense packings 
in the Euclidean spaces. Via the canonical Q-embedding of arbitrary number field K into 
both the prime ideal p and its residue field k can be embedded as discrete subsets 
in Thus we can concatenate the embedding image of the Cartesian product of n 

copies of p together with the image of a length n code over k . This concatenation leads to 
a packing in the Euclidean space R”[^'®. Moreover, we extend the single concatenation to 
multiple concatenation to obtain dense packings and asymptotically good packing families. 

Eor instance, with the help of Magma, we construct a 256-dimensional packing denser than 
the Barnes-Wall lattice BW 256 . 
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1 Introduction 

The classical problem of packing non-overlapping equal spheres densely in an n-dimensional Euclidean space 
has attracted the interest of numerous mathematicians for centuries. Many methods and results from different 
disciplines, such as discrete geometry, combinatorics, number theory and coding theory, etc. have been 
involved in this problem, while some explicit fascinating dense constructions and asymptotically good packing 
families have been found. For a detailed survey on the development in this field, the reader may refer to the 
books of Cassels [B], Conway and Sloane [5], Zong [50]. If the centers of the packed spheres form a discrete 
additive subgroup of R" (lattice), we call it a lattice packing. 

Sphere packings are continuous analogues of error-correcting codes in the Hamming space |0]. The basic 
problem of error-correcting codes is to seek the maximum size of a code given length, alphabet size and 
minimal Hamming distance, in other words, dense packing of points in Hamming space such that each pair 
of distinct points is separated at least by the minimum Hamming distance. In the digital world, error- 
correcting codes are widely employed in information storage and transmission, for example, the blue ray 
storage format and the communication between space stations and the Earth. Based on the similarities 
between sphere packings and error-correcting codes, the results in sphere packings can potentially contribute 
to the development of error-correcting codes. 

Meanwhile, some constructions of dense lattice or non-lattice packings are inspired by constructions in 
coding theory. For example, similar to concatenated codes, Leech and Sloane’s “Construction A” method 
concatenated certain binary codes together with 2 • Z” to construct new lattice packings (see details in 
[5] Chapter 5]). Bachoc |5] generalized the method to construct modular lattices using codes over finite 
involution algebras. In Construction C [201 Chapter 5], the binary expansion of the coordinates in Z" is 
considered. A point is a packing center if and only if its first £ coordinate arrays are codewords in certain i 
binary codes respectively. 

Let uj = [IS] further investigated the concatenating method. Instead of packings in Z", he 

considered the packings in where Ok = ^ [w] denotes the ring of integers in the number field Q(\/^). 
That is, for a non-zero prime ideal of Ok, fP” can be embedded as a lattice L in R^" via the canonical 

^This research began when the author was a PhD candidate at Nanyang Technological University. 
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Q-embedding of the special number field Q(-\/^) into Hence any subset V of fp" can be regarded 
as a packing in R.^". Then he replaced the binary expansion in Construction C by the fp-adic expansion, 
and concatenated the lattice L with some special codes over Ok/'^- This method produces several dense 
packings in low dimensions attaining the best-known densities and also obtains an unconditional bound on 
the asymptotic density exponent A ^ —1.265 (see m)- Cheng [7] applied this concatenation to multiplicative 
lattices and improved the asymptotic density of packing families derived from multiplicative lattices. One 
natural question is whether we can use arbitrary number field instead of the special quadratic number field 
Q(-v/^) to generalize the constructing method. 

In this paper, we extend Xing’s method to a general level, i.e. we employ the ideals in Ok instead of 
Z[a;], where Ok denotes the ring of integers in an arbitrary number field K. Suppose the extension degree 
of K over Q is m. Minkowski interpreted the elements in K as points in the m-dimensional Euclidean space 
R™. The interpretation is called “Minkowski Theory” in algebraic number theory (see [H Section 1.5] or 
[l8l Section 5.3]). 

By Minkowski’s interpretation, we can use the canonical Q-embedding to construct lattices in R"* from 
ideals in Ok- In this way, for a non-zero prime ideal p of Ok, the Cartesian product of n copies of p can be 
embedded as a lattice in R™". Meanwhile, the codes defined over OkIP with length n can also be embedded 
as finite subsets of R™”. Then we can proceed with the concatenating method as Construction A on these 
two kinds of subsets of R™" to construct dense packings. For instance, we construct a 256-dimensional 
packing with center density 6 satisfying log 2 6 ^ 208.09, which is larger than 192 of the Barnes-Wall lattice 
BW 256 (see the table of dense packings in jS] Table 1.3]). Furthermore, for different choices of number helds 
and prime ideals, we also provide several asymptotically good packing families. 

In sectionjT] we recall the formal definitions and necessary properties of sphere packing densities, algebraic 
number fields and codes. They play important roles in the main results. In section [3] we describe our 
generalized concatenating method in detail. Several examples of dense packings and asymptotically good 
packing families for different choices of number fields and prime ideals are provided. The detailed numerical 
results are provided in Tables [T] - |4j In Section H] we conclude our main contribution. 


2 Preliminaries 


2.1 Sphere Packing 

Let V be the set of centers of packed spheres and Bm{R) be the set 


|v = (oi, • • • jOat) e R^ : 





where ||v|| denotes the Euclidean length of vector v. As a sphere packing construction is determined by the 
arrangement of the sphere centers, we just use V to denote the corresponding packing. 

For a packing V, the radius of the equal packed spheres is r = d£:(P)/2, where dE{V) is the minimum 
Euclidean distance between two distinct points in V. Then the density A(7^) of packing V is defined as 


A(7^) = limsup 

R—^oo 

where Vm is the volume of the unit sphere in R’ 

r 


\VC^BN{R)\ 


yo\{Bn{R 
that is 


Vn 


(A/2)!’ 

2^7r(^-i)/2((A_l)/2)! 

m 


■r^ -Vn 

+ 0 ) 


if N is even; 
if N is odd. 
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The sphere packing problem is to construct packings obtaining large density A(P). Moreover, the center 
density SiJ’) and density exponent A(P) are defined respectively as 

<5(7’) = ^, A(P) = llog2A(P). 

liV = L forms a lattice, the density of lattice packing L can be simplified as 

{dE{L)/2)^VN 


A(T) = 


det(L) 


where det(L) is the determinant of L. Note that from Minkowski’s Convex Body Theorem (see [131 Theorem 
1.4] or [lH Theorem 4]), for any rank N lattice L, the minimum Euclidean distance of L satisfies 


dE{L) < VAdet(L)i/^. 


( 1 ) 


When we explore the asymptotic behavior of a packing family T = as the dimension N of packing 

tends to oo, we consider the asymptotic density exponent of the family 


A(J-)=limsup4log2A(iPW). 

iV V / 


( 2 ) 


Note that by Stirling formula, as —>■ oo, we have 


N N 1 

log 2 Vn = -y log 2 ^ “ 2 - e, 


(3) 


where 0 < e < (log 2 e)/(61V). 

For the asymptotic aspect, Minkowski gave a nonconstructive bound that asserts that there exists some 
packing family J" such that the asymptotic density exponent X{J-) satisfies A(J^) > —1 (see [6l p.l84]). It 
is of interest to construct families with A(J^) > —oo explicitly (such families are called asymptotically good). 
Known constructive bounds for families with polynomial or exponential construction complexity in terms of 
dimension N are listed in the book of Litsyn and Tsfasman [T51 p.628]. Rush [T7] proved that Minkowski’s 
bound on asymptotic density A > — 1 can be attained by lattice packings constructed from codes. More 
recently, Gaborit and Zemor m improved the density by a linear factor to the quantity of the form cn2 " 
for constant c. 


2.2 Algebraic Number Fields 

Let K be an algebraic number field of degree m over Q, and let Ok be its ring of integers. Suppose K has 
s real embeddings 

pi, • • • ,ps : K 

and t pairs of complex conjugate embeddings 


CTi,cri,--- K ^ C. 

Thus m = s + 2t. We consider the canonical embedding t : K ^ that is, for any a G K, r(a) is the 

following vector 

(pi(a),-- - , Ps{oi),'/2iSiai{a),\/2^ai{a),-■ ■ ,'/2ifiat{a),\/2^at{a)'^ , (4) 

where 5R denotes the real part and 3 denotes the imaginary part of a complex number, r can be directly 
extended to AT" M™", which is also denoted by r in this paper without causing any confusion. 

Note that in the software Magma V2.21-4 00, the embedding r defined above is called the Minkowski 
map and it is employed to compute the “minimum norm” (square of minimum Euclidean distance) of a 
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Minkowski lattice. In Subsection 13.41 we will use the software to compute the density of some explicit 
examples from our construction. 

Set Tr = Trace^f/Q and N = Norm^^/Q. By the Minkowski Theory [TH Section 1.5], the embedding 
r maps the non-zero ideals of Ok to some special lattices in R."*. The following lemma characterizes the 
determinant of such lattices. 

Lemma 2.1 ( \1A[ Section 1.5]). For any non-zero ideal I C Ok, Lj := t{I) is a lattice of rank m. The 
determinant of Lj is 

det Li = Af{I)y/\D^\, 

where Dk is the discriminant of K andAf{I) = [Ok ■ I] is the absolute norm of I. 

Generally it remains hard to compute the minimum Euclidean distance of lattices. However, we can 
estimate a lower bound on the Euclidean length of non-zero points in the lattice Lj. The proof is sketched 
here and a similar discussion on general ideal lattices can be found in [3]. Note that in [5], the distance 
between two points is defined using some special quadratic forms, while in this paper, we only focus on the 
standard Euclidean distance, which refers to the original meaning of packing in the Euclidean space and 
brings convenience in computation (Magma uses standard Euclidean distance as built-in measure). 

Lemma 2.2. Let I C Ok be a non-zero ideal. For any 0 a S /, the Euclidean length of the vector 
T{a) G Li satisfies 


||T(a)|| ^ y/m- |N(a)|^'^'" ^ ^/m- 

In other words, dsiLi) ^ y/rri ■ . 

Proof. The Euclidean length of the vector T{a) satisfies 

S t 

lk(a)f = + 




^ (s -|- 2t) 


■ Y[(Tj{a)aj{a) ■ W(Tj{a)aj{a) 

l/n 


= m (N(a))" 

The last ^ becomes an equality if and only if the ideal / is a principal ideal generated by a. 


□ 


The following lemma is a basic fact in algebraic number theory (see (M] Section 1.2]). 

Lemma 2.3. For any non-zero element (3 G Ok, we have |N(/3)| ^ 1. 

From the above lemmas, the minimum Euclidean length of non-zero elements in Lok satisfies ^^(Lok) ^ 
y/rri. Indeed, as the vector r(l) has Euclidean length y/rn, we have 


dsiLoK) = Vm. 


( 5 ) 
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2.3 Coding Theory 

We recall some notations and results in coding theory. 

For a g-ary code C, let n{C),M{C) and dniC) denote the length, the size, and the minimum Hamming 
distance of C, respectively. Such a code is usually referred to as an (n(C'), M(C'), d//(C'))-code. Moreover, 
the relative minimum distance g{C) and the rate R{C) are defined respectively as 


g{C) = 


dniC) 

n{C) ’ 


R{C) = 


log,M(C) 

n{C) 


In particular, if a code C forms a linear space over Fg, then the code C is called an [n(C), fc(C'), d_f/(C')]- 
linear code over Fg, where k{C) := log,j M(C') is called the dimension of C. In this case, the rate R{C) = 
k{C) 
n{cy 

Let Uq be the set of the ordered pair {g,R) S for which there exists a family of q-ary codes 

such that n(Ci) increasingly goes to oo as i tends to oo and 


g = lim p(C'i), and R — lim R(Ci). 

2—)-00 i—¥00 


Here is a result on Uq-. 


Proposition 2.4 ([T51 Section 1.3.1] or m Proposition 3.1]). There exists a continuous function Rq{g), 
g € [0,1], such that 

Uq = {{g,R)GM.^ lOifRifRqig), 0 < p ^ l} . 

Moreover, Rqi^) = 1, Rqig) = 0 for g € [{q — I)/?, 1]; and Rq(g) decreases on the interval [0, {q — f)/q\. 

For 0 < p < 1, the q-aiy entropy function is given as 


Hq{g) = glogqiq - 1) - plog,j p- (1 - p)log,j(l - p). 

The asymptotic Gilbert-Varshamov (GV) bound indicates that 

Rq{g) ^ Raviq, g) ^ - Hq{g), for all p e ^0, . (6) 

Moreover, for any given rate R, there exists a family of linear codes which meets the GV bound (see m 
Section 17.7]). 


3 Main Results 

In this section, we introduce our new construction of dense sphere packings. In particular, several explicit 
constructions and asymptotically good packing families are provided at the end of this section. Our idea is 
to apply special concatenating methods on Lj defined in Lemma l2.11 

Let [if : Q] = TO and let p be a prime ideal of Ok- Assume that the residue field Fp = Ok/P of p is 
isomorphic to the finite field Fg. Let r be the canonical embedding defined in Eq. O- 

For i S N, let Uy denote the Gartesian product of n identical copies of lattice Lpi = r(p*). Then L” is 
a lattice of rank mn by Lemma l2.II Moreover, the determinant satisfies 

detuy = (detLpi)", 

which is straightforward from the definition of lattice determinant. In addition, by the definition of L^, the 
minimum Euclidean distance satisfies 

dEiL;y = dE{Lpi). 
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For simplicity, in this section we write Li = Lpi for short. 

As Ok is a Dedekind domain (see [Ml Section 1.3] or [T1 Chapter 9]), let k = OkIP, then we have 
dim„pYp*+^ = 1 for all i G N. We know k = F^. Thus for i G N, we can choose the sets Si := 

= 0, , • • • , Qfq ^ I C p* such that 

af mod p*+\ • • • , af mod p^+^ (7) 

represent the q distinct elements in pYp*"*"^. 

3.1 Concatenation with One Code 

We hx one index io G N^i in the discussion here and generalize the result to a family of indices in the next 
subsection. Let 7^ be a subset of which can be regarded as a packing in R"*". 

Given a g-ary code C with length n, in order to perform our concatenation, we take the code alphabet set 
of C to be Sig-i C p*o“i (see Eq. (I?])). In this way, C can be regarded as a hnite subset of O^. The advantage 
of using 5^0-1 as the alphabet set is that it can help us bound the packing radius of our construction. The 
details are included in the proof of Proposition 13.21 

We define r(c) for each codeword c = (ci, • • • , c„) G C C as 

r(c) := (t(ci), • • • , r(c„)) G R"*", 

and take t{C) := {t(c) : c G C}. We consider the concatenation 

t{C) + P := {a + p G R™" : a G r(C), p G P} , 

which corresponds to a packing in R"*”. We will analyze the density of such a packing in this subsection. 

The following lemma characterizes the minimum Euclidean distance of the concatenation r(C) + P, which 
will play an important role in later discussions. 

Lemma 3.1. Let p be a non-zero prime ideal in Ok with absolute norm A/'(p) = q {i.e., the residue field 
Fp = Fg). For zq G let 

(i) V be a subset of Lf^, 

(a) C be a q-ary (n, M,dc)-code over the alphabet set ^,0-1- In addition, C contains the zero codeword. 
Then the minimum Euclidean distance of the packing t{C) + V satisfies 

ds (t(C') + V)fi min|dE(Lio_i)v^, d£;(7’)| . 

In particular, if dE{Li^-i)^/dc ^ de {V), we have exactly 

dE {t{C) +V) = dE{V). 

Proof. Let t(c) + p be a non-zero element in t(C') -I- V, where c G C and p G 7^. If c = 0, then 

lk(c) +p|| = IIpII ^ dE{P). 

If c 7 ^ 0, without loss of generality, we assume c = (ci, • • • , Ce, 0), where e is the Hamming weight of c 
and Ci 7 ^ 0 for 1 ^ z ^ e. As V C L", we can further assume p = (r(pi),--- ,t(p„)), where pj G p®° for 
1 ^ j ^ Thus 

t(c) + P = (r(ci -k pi), ■■■ , T{Ce + Pe),T{pe+l), ■■■ , t(p„)) . 
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As Cl, • • • , Ce are non-zero elements in Sig-i, i.e., for 1 ^ i ^ e, ci S ^ \ P*°j we have that Ci + pi are 
non-zero elements in p®°“^ for 1 ^ ^ e. Thus 

e 

\\t{c) + pf ^ ^ ® ■ d%{L^g-l) ^ dc ■ d|(Lio-i)- 

i^l 

In summary, the minimum Euclidean distance of the packing t{C) + V satisfies 

dsiTiC) +V) ^ mm^dE{L^g-i)-/d^, dE(P)| . (8) 

In particular, if dE{Lig-i)^/dc ^ dE^V), then from Eq. ([8]) above, we have the minimum distance 
dE {t{C) -fP) > dEiJ’)- Meanwhile, as C contains the zero codeword, P is a subset of t{C) -\- P, thus 
dEiP) ^ dE (t(C') -|- V). Finally we have exactly 

dE {T{C)+V)=dE{P). 


□ 

The following proposition generalizes Proposition 2.3 of m, where only the special case K = Q(\/^) 
is discussed. 

Proposition 3.2. Under the same assumption on P and C as in Lem,m,a \S.l\. 

(i) if dE{Lig-i)\/dc ^ dE iP), then the density of t{C) +P as a packing in R™” is at least M ■ A (P), 
where A {P) is the density of the packing P as a packing in M™”, and M is the size of the code C; 

(a) ifP is a lattice and C satisfies that for any u, 0 S C, the sum t(u) -|- T(ti) is equal to T(rD) -I- p for some 
to G C and p € P, then t(C) +P is also a lattice. 

Proof, (i) We denote the volume of the unit sphere in by Vat and the sphere of radius 6 by ,6 at (6). Let 
s = max{||r(c)|| : c S C}. For any c € C and p € P D Bmn{h), we have t(c) -f p G (t(c) -I- P) fl Bmn{b + s). 
This implies that 


|(t(c) PP)UBmn{h + s)\ ^ \PUBmu{b)\ ■ 

Moreover, as the elements in the alphabet set Sig-i of C represent the q distinct elements in p'o^^/p'o and 
P U L” , we have if Ci ^ Cj then 


(r(c,) +P)n (t(cj) -f P) = 0. 


We write d = dE {t{C) -I- P) for short and from Lemma [3.II we have d = dE{P). Finally we obtain 




A(r(C)+P) 

\{T{C)+P)UBmu{b + s)\{d/2r'^Vmu 
YO\{Bmuib + S + d/2) 


lim sup 
6^00 

lim sup 

b—^oo 


(Ecec 1(^(0 + ^) n Bmn {b + S)|) {d/2r^ K, 

vo\{Bmn{b + S + d/2) 

|C| • IP n (d£(P)/2r” 


vol {Bmn (b - 
\C\- A{P) = M ■ A{P). 


■d/2)) 


(ii) By the definition of a lattice, this part is true. 


□ 
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3.2 Concatenation with a Family of Codes 

For a family of g-ary codes where £ € N^i, we can take the code alphabet set of C to be Si (see 

Eq. ([7])). Note that the choice of the alphabet set of the codes is only used in the proof. For the computation, 
we only care about the length, size and minimum Hamming distance of the codes. 

Similar to Subsection ixn for any subset V of L^, and a family of g-ary codes C = length 

n, we consider the concatenation 

t-i (i-i 'j 

r(C)+iP:=^r(C)+iP= + :a, Gr(C),pGiPL (9) 

i=0 I i=0 J 

Note that for each 1, as the code alphabet set of Ct is Si, the set r(Ci) is a finite subset of L". 

Meanwhile 7^ is a subset of L". As for any 0 ^ ^ t', the lattice L" is a sublattice of the concatenation 

Eq. ([5|) makes sense within . 

Proposition 3.3. Let p 6e a non-zero prime ideal in Ok with absolute norm M{p) = q. For £ G Nj.i, let 


fi) 

(ii) 


V be a subset of L'l, 


C = ® family 

r4(P) 


and dci ^ 


dULi) 


In 


of q-ary codes, where Ci is an {n, Mi,dci)-code, the alphabet set of Ci is Si, 
addition, for each 0 ^ i ^ £ — 1, Ci contains the zero codeword. 


Then the density of t{C) -\-V as a packing in R™" is at least A(P) • Y[i=o 
Proof. Let Vg :='P C L". For i from £ — 1 to 0, we can recursively define 


V, :=r(C)+Pm CL-. 


Note that Vq = t(C) + V. 

We claim that for 0 ^ f ^ — 1, 


dsiPi) = dEiVi+i) and A{V,) ^ Mi ■ A{'P,+i). (10) 

We use induction on fc = £ — i to prove the claim. When fc = l,f = £—1, as the minimum Hamming 
distance of Ce-i satisfies dc,,_^d\[Li-i) ^ d\{fP) = d\{fPi), by Lemma [5?T] and Proposition 13.21 fil. we have 

dE(fPe-i) = dE (r(C^_i) + Vi) = dEiVi) and A(P^_i) > • A{'Pt). 


Suppose for all k in the range l^fc<fco + l^£, the statement Eq. (nni is true, i.e., i = £ — k, 

dEiVt-k) = dEiVt-k+i) and A{Ve-k) Mg_k ■ A{Ve-k+i) (11) 

is true. By induction, we need to prove 

dEiVt-ko-i) = dE{Pe-ko) and A{fPt-ko-i) ^ ' A{'Pt-ko)- 

From Eq. m, we have 

dEiVi-ko) = dE{Fi-ko-i-i) = ■ ■ ■ = dEiVi-i) = dE{Vi) = dEiV). 

Besides, the minimum Hamming distance of Ct-ko-i satisfies 


dci-ka-id\{Li-ko-i) ^ d%{V) — d\{Vi.-ko)- 





Then by Lemma ET] and ProDOsition l3.2l fib as Vt-ko-i = T(Q-feo-i) +^^-fco) get 


dEi'Pe-ko-i) = dEi'Pi-ko) and A{Vi-ko-i) ^ Mi^ko-i ' ^{'Pt-ko)- 

Thus we have proved the claim Eq. (fTUl) . From the claim, we easily obtain A(Vo) ^ ^(P) ■ Ilfco ^ 

Remark 3.4. From the above proposition, given a dense packing V from the canonical embedding Eq. o 
on some prime ideal in a certain algebraic number field, we can concatenate several codes satisfying the 
requirements in Provo sition \S.t^ to V to obtain some new denser packings. The density increases by a ratio 
larger than the product of the sizes of these codes. 

Note that generally it is hard to determine the minimum Euclidean distance of P C C R™". However, 
we can consider special choices of P with some algebraic structure. The following corollary, which plays a 
crucial role in Subsection 13.41 for Examples 13.81 - 13.111 considers the case P = L'f. The advantage is that 
for all i £ N, we have dE{Lf) = dE{Li)- Instead of searching for the minimum Euclidean distance of Lf in 
we can search for the minimum Euclidean distance of Li in R™. In our examples, m is small such that 
Magma can be used to compute the minimum Euclidean distance. 


Corollary 3.5. Let p be a non-zero prime ideal in Ok with absolute norm M{p) = q. For i £ N^i, let 
C = {Ci}i^Q be a family of q-ary codes, where Ct is an (n, Mi,dCi)-code, the alphabet set of Ci is Si, and 

.In addition, for each 0 ^ i ^ — 1, the code Ci contains the zero codeword. 


dCi > 


dULi) 


Then the density of t{C) + L'f as a packing in R™" satisfies 

(dE(L^)/2)™"' Vmn ^ ^ 


A{r{C)+L'l)^ 




2=0 


where Dk is the discriminant of K and Vmn is the volume of the unit sphere in Moreover, the center 

density S = S{r{C) + L^) satisfies 

i-i 

log 2 5 ^ mnlog 2 dEiL^) -mn- n^\og 2 q - - log 2 \Dk\ + ^ log 2 


2=0 


In particular, if for all 0 ^ i ^ £ — 1, the code Ci is a linear code of dimension ki = log Mi, then 


e-i 


log2 S ^ mnlog2 dE{Li) - mn - nilog2 9 ^ log2 \^k\ + log2 q-^h. 


2=0 


Proof. Note that M{p^) = q^ and then by Lemma I^TT] 


detL^ = (detL^)” = {q^ 


□ 


Remark 3.6. For concatenation with finite codes, we care more about the dimension of the packing as we 
need to compare the new construction with the known good packings in the Euclidean space with the same 
dimension. Hence we fix the length n of the codes first. As for 0 ^ i ^ £ — 1, the lattice Li is a sublattice of 

Lq, we have dE{Li) ^ dE^Lo). Thus the length n satisfies n ^ ^ by Eq. ([5|) and Lemma 

?. This bound will be used in the 


Therefore, we can concatenate at most £ = — log n 


computation of Subsection \3.4 


L 2 


4M' 

d^iLo) 
codes to L 
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3.3 Asymptotic Properties 

In this subsection, we show that our construction will lead to asymptotically good packing families, which 
means we can obtain several constructive bounds on the asymptotic density exponent Eq. ©■ 

We will employ families of linear codes that meet the GV bound. As there exist only exponential time 
algorithms or randomized polynomial algorithms to construct such families (see [T^ Section 17.7]), our 
corresponding asymptotic density exponent bounds are exponential constructive bounds. 

From Eq. o and Lemma 12.21 for 0 ^ i ^ — 1, we know 




( 12 ) 


Different from the finite concatenation in Subsection l3.21 for asymptotic results, we first fix 
the packing for each £ G N>i. The length ni of the codes is set depending on £. 

i|( 


We uniformly set ng = . Then by Eq. (IT^ . ng ^ 

on the GV bound Eq. for 0 ^ z ^ — 1 , we can choose g-ary 


dl{Lg) 


dl{L, 


for all 0 ^ z 


£ and construct 
^ £ — 1. Based 



■4MT 

d%{Li) _ 


linear code 


such that the rate is 


ng 


^Rcviq, = 


where the relative minimum distance satisfies 


(ne) _ 
i 


' dliLgy 

d%{Li) 

ng 


Let 


ein,) |(7 




z . For each ng, define a packing 
J i=0 


V 


("z) . 


:= r(C("^)) + L; 


as in Eq. (0. 

Next we take £ increasingly to oo. The following proposition describes the asymptotic density exponent 
of the packing family where ng = tends to oo as £ goes to oo. 

Proposition 3.7. The asymptotic density exponent of the family T satisfies 

1 1 771 


+ limsup ( log 2 dEiLg) - i log 2 ng - 

Z-i-no \ 2 77Z 




("z)\ 


(13) 


2 = 0 


where H'^(g) = Hq(g) for Q < g < - - and Hq{g) = 1 for - - ^ p ^ 1. 
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Proof. By Corollary 13.51 


By Eq. ([3]), 
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({dE{LdirV. 


i-i 


\{T) ^ limsup-log 9 , , ^ , 

mm y 

= -1 - :r^l0g2 \Dk\ + 


mng_ n 

i=0 


2m 


1 


+ limsup log 2 dE{Ll) -log 2 q H-log 2 Vr 


l—¥00 


limsup 

—^OO \ 


■ l0g2 Vr 


2 * mni 


.. / 11 
limsup --log 2 


E 

i=0 

1 


mni 


^ log 2 q • 


E 


l0g2<?-fc^)' 


mui 




27re 2mni 


mn£ I 

log 2 q ^ 
log2 mneTT H-—— y .' 


1 , m / 1 , log 2 q 

^ -olog 2 :^-hlimsup --log 2 n^H- 

z Z7re ^-s-oo \ ^ w 


“o ^og2 -^ limsup -- log2 Hi H- 

2 27re e^oo \ 2 m 


2 27re i 
In summary, we get Eq. (USD. 


1 . , r Ml 1 1 log2 9 

“o log 2 -^ limsup — log 2 q - - log 2 ne - 

\ m 2 m 


m Tie / 

i=0 / 

'^Rgv ( 9 , £>1""^) 

z— 0 X 

i: (1 - "i (ft 


i=0 


□ 


3.4 Examples 

We use some explicit examples to illustrate our new construction. Here Magma V2.21-4 011] will be employed 
to obtain the numerical results. 

For the finite concatenation, we use the linear codes from m which have an explicit construction. For 
the asymptotic result, we set sufficiently large £ to approximate the limit. 

Example 3.8. Consider the number field K = Q[a], where a is a root of the irreducible polynomial 

f{x) = x'^ — + X + 1 £ Q[x]. 

iiT/Q is a totally complex number field with extension degree [if : Q] = 4 and absolute discriminant \Dk\ = 
117. The absolute discriminant is the smallest one of all totally complex number fields with degree 4 (see 

m)- 

(i) We consider a prime ideal p lying over 3 S Z. The Magma code of this example is listed in Appendix 
HI while for other examples, the Magma code can be modified from this template. 

The key outputs of AvvendixlA] are listed here: 


The degree of K=Q[x]/(x~4 - x"3 - x~2 + x + 1) is m=4; 
The absolute value of the discriminant of K is |d|=117; 
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P is a Prime Ideal of 0 
Two element generators; 

[3, 0, 0, 0] 

[2, 1, 1, 0] lying over 3 with norm q=9; 


The above statement means p = (3, 2 + a + a^). 

Finite Concatenation: 

We can concatenate at most 3 linear codes to the lattice 
constructed by 64 copies of L_{P~3}, 
whose Hamming weights are required respectively at least 
27;9;3; 


Referring to m, as the norm of p is 9, we can use the existing 9-ary linear codes Co,Ci,C 2 with 
parameters 

[64,25,27], [64,49,9], [64,61,3] 

respectively. The sum of the dimensions is 135. Considering the concatenation in Corolla,ru 1,9.5L we 
obtain a packing with dimension 4 x 64 = 256, whose center density S satisfies 

Our packing is in dimension 256 with Log_2(center density) 
at least 208.088204168043224246772217517 


Note that our packing is denser than the Barnes-Wall lattice BW 25 e, whose density is listed in the table 
of dense sphere packings la Table 1.3]. 

Using this ideal to construct a packing family as in Proposition \3. 7[ we get the following result. 

Asymptotic result: 

The asymptotic density exponent of the packing family is 
at least -1.442426720 


Note that the above result means our packing family is asymptotically good. 

(ii) For this field, we similarly analyze other prime ideals and list part of our numerical results in Table\^ 
where S is the center density and A is the asymptotic density of the corresponding packing family. 

Example 3.9. Consider the number field K = Q[q!], where a is a root of the irreducible polynomial 

f{x) = — 2x — 1 G Q[x]. 

K/Q is a totally real number field with extension degree [iG : Q] = 3 and absolute discriminant \Dx\ = 49. 
The absolute discriminant is the smallest one of all totally real cubic number fields (see m)- The numerical 
results on our examples are partially listed in Table\^ 
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Table 1: Examples constructed from Q[a;]/ — a;^+a;+l) 


p 

9 

dimension 

log 2 6 at least 

A at least 



180 

108.52 


(3, 2 “h Of -|- Of^) 

9 

256 

208 .091 

-1.442 



512 

590.52 


(7,4 +a) 


256 

190.63 


7 

400 

410.15 

-1.453 


^ The packing can be explicitly constructed with log 2 5 larger than 192 
given by the Barnes-Wall lattice BW256- 


Table 2: Examples constructed from Q[a:]/ (a:^ + — 2a: — l) 


P 

9 

dimension 

log 2 S at least 

A at least 

(2) 

8 

255 

134.46 

-1.628 

(7, 5 +a) 

7 

192 

255 

83.68 

157.63 

-1.585 


Example 3.10. Consider the number field K = Q[a], where a is a root of the irreducible polynomial 

f{x) = — 1 G Q[x]. 

K/Q is a number field with extension degree [E' : Q] = 3 and absolute discriminant \Dk\ = 23. The absolute 
discriminant is the smallest one of all cubic number fields (see m)- The numerical results on our examples 
are partially listed in Table\^ 


Table 3: Examples constructed from Q[x]/ (x^ + x^ — l) 


P 

9 

dimension 

log 2 6 at least 

A at least 

(2) 

8 

96 

192 

24.70 

115.40 

-1.429 

(5, 2 +a) 

5 

150 

180 

69.47 

101.01 

-1.445 

(7,11 +a) 

7 

255 

187.32 

-1.430 


Example 3.11. Consider the number field K = Q[a], where a is a root of the irreducible polynomial 

/(x) = X® + x^ + 1 G Q[x]. 

K/Q is a number field with extension degree [ET : Q] = 6 and absolute discriminant \Dk\ = 19683. The 
polynomial /(x) = x® + x^ + 1 is the 9 th cyclotomic polynomial over Q. The numerical results on our 
examples are partially listed in Table(^ 

4 Conclusion 

This paper provides a new method to construct dense packings using canonical Q-embedding of algebraic 
number fields, and special codes over the residue field of some prime ideals. With the help of Magma V2.21-4, 
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Table 4: Examples constructed from Q[a;]/ {x^ + + l) 


p 

9 

dimension 

log 2 5 at least 

A at least 

(3, 2 +a) 

3 

180 

192 

109.71 

122.72 

-1.868 


several explicit constructions were provided in Examples 13.81 - 13.111 Especially, in a packing denser 

than the Barnes-Wall lattice BW 256 was obtained. Moreover, this method can be utilized to construct 
asymptotically good packing families. For each number field and prime ideal in Tables [T] - |4l a lower bound 
on the asymptotic density exponent of the corresponding packing family was provided. 
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A Magma Code for Example 13.8 

R:=RealField(10); Q:=RationalField(); 

W<x>:=PolynomialRing(Q); f:=x“4“x“3-x“2+x+l; 

/*Input a polynomial over rational number field Q*/ 
if Isirreducible(f) then /* Test whether f is irreducible*/ 

K<a>:=NumberField([f] ); 0:=MaximalOrder(K); m:=AbsoluteDegree(K); 
printf'The degree of K=Q[x]/(y,o) is m=y,o; \n" ,f ,m; 
d:=AbsoluteDiscrimincLnt(K); 

printf'The absolute value of the discriminant of K \ 
is I d I =y,o; \n" ,d; 

p:=3; /*Choose p=3 as a base prime number*/ 

J:=Decomposition(0,p); P:=J[1,1]; q:=Norm(P); 
printf'P is a y«o lying over y,o with norm q=y,o;\n" ,P,p,q; 

printf"\nFinite Concatenation:\n"; 
n:=64; 1:=Floor(m/2*Log(n)/Log(q)); 

/*Set the length of the codes*/ 

L: =MinkowskiLattice(P'‘l); b:=Minimum(L); 
printf "We can concatenate at most y,o linear codes to \ 
the lattice constructed by y,o copies of L_{P'‘y,o}, \n \ 
whose Hamming weights are required respectively at least \ 
\n",l,n,l: 

for t in [0..1-1] do 

Z: =MinkowskiLattice(P'‘t) ; h:=Ceiling(Minimum(L)/Minimum(Z)) ; 
printf "y,o;" ,h; 
end for; 
printf"\n"; 
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printf"Refer to the Codetable.de to get the largest dimension \ 
of the corresponding linear codes with length 7,o; \n",n; 
fromtable:=135; 

printf "Here the sum of dimensions T=7«o;\n" ,f romtable; 

/*For this example, the sum of the dimensions is 135*/ 
c : = (m*n*Log(b'‘0.5) -n*l*Log(q)-0.5*n*Log(d) )/Log(2) -m*n; 
printf'Our packing is in dimension y»o with \ 

Log_2(center density) at least\n 7,o\n",n*m,\ 
c+Log(q)/Log(2)*fromtable; 

printf"\nAsymptotic result:\n"; 

1:=1000; n:=Ceiling(q“(2*l/m)*d‘(l/m)); 

/*Set 1 sufficiently large to approximate the limit*/ 

L: =MinkowskiLattice(P'‘l) ;b: =Minimum(L); 

Sum:=0; 

for t in [1-1..0 by -1] do 
Z: =MinkowskiLattice (P'‘t); 
g:=Ceiling(Minimum(L)/Minimum(Z)); 
varrho:=g/n; 

if varrho le (q-l)/q then 

Sum:=Sum+(varrho*Log(q-1)-varrho*Log(varrho)\ 

-(1-varrho)*Log(l-varrho))/Log(q); 
else Sum:=Sum+l; 
end if; 
end for; 

Lambda:=-1-1/(2*m)*Log(d)/Log(2)-0.5*Log(m/(2*Pi(R)*Exp(l))) \ 
/Log(2)+Log(b‘0.5)/Log(2)-0.5*Log(n)/Log(2) \ 

-l/m*Log(q)/Log(2)*Sum; 

printf'The asymptotic density exponent of the packing family is \ 
\n at least y»o",Lambda; 

else 

printf'The polynomial ‘/o is not irreducible over Q.\n",f; 
end if; 
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